Mon. Not. R. Astron. Soc. OOP. [TlfTil COOOO') 



Printed 4 February 2008 



(MN I*TeX style file v2.2) 



Discovery of an open cluster with a possible physical 
association with a planetary nebula 



. C. Bonatto\ E. Bica^ and J.F.C. Santos Jr.^ 

, ^ Departamento de Astronomia, Universidade Federal do Rio Grande do Sul, Av. Bento Gorifalves 9500 
• Porto Alegre 91501-970, RS, Brazil 

' ^ Departamento de Fisica, ICEx, Universidade Federal de Minas Gerais, Av. Antonio Carlos 6627 
' Belo Horizonte 30123-970, MG, Brazil 



4 February 2008 



ABSTRACT 

We report the discovery of a new open cluster (OC) in the Galaxy at £ — 167.0° 
and b = -1.0°. Its field includes the planetary nebula (PN) PK 167-0.1. We study 
the possible associations of the PN/OC pairs NGC 2818/NGC 2818A, NGC 2438/M46 
(NGC 2437), PK 6+2.5/NGC 6469, as well as of the PN PK 167-0.1 with New Cluster 1. 
The analyses are based on near- infrared colour-magnitude diagrams (CMDs) and stel- 
lar radial density profiles (RDPs). NGC 6469 is located in a heavily contaminated 
bulge field. The CMD morphology, especially for the latter two cases, is defined with 
a field star decontamination algorithm applied to the 2MASS J, H, and Kg photome- 
try. Field decontamination for the DCs NGC 2818A and M46 produced better defined 
CMDs and more accurate cluster parameters than in the literature. Those pieces of 
evidence point to M46 as physically associated with the PN NGC 2438. The same oc- 
curs for the OC NGC2818A and the PN NGC 2818, however previous radial velocity 
arguments indicate that they are not associated. The OC NGC 6469 does not appear 
to be associated with the PN PK 6-1-2.5, which probably belongs to the bulge. Finally, 
the distance of the OC New Cluster 1 is consistent with a physical association with 
the PN PK 167-0.1. 

Key words: (Galaxy:) open clusters and associations: general; (ISM:) planetary 
nebulae: general. 



1 INTRODUCTION 

Planetary nebulae (PNe) are late stages in stellar evo- 
lution, occurring for star s with mass in the range 1 — 
6.5 M(7> (| WeidemannI |2000| '). possibly with the upper limit 
at ^ SMq. Most of the difficulty in understanding such 
stages that lead to PNe is associated with the determina- 
tion of the PN distance and its progenitor mass. Probably, 
the largest source of error is the distance which, 
in m ost cases, is known with a ~ 50% uncertainty 
(e.g. IZhand or more. However, both prob- 

lems can be minimised if the PN is physically as- 
sociated with a star cluster. Star cluster distances, 
in general, can be determined with a high preci- 
sion and, in the case of physical association, the 
PN progenitor mass can be assumed to be ~ 10% 
large r than that of the turnoff (TO), stars referred 
to by iMevnet. Mermilliod &: Maederj (| 19931 ) as red- 
turnoff objects. Exceptions are the cases where mass 
transfer from a close binary is involved. In fact, 
recent radial velocity measurements, although re- 



stricted to relatively small samples, suggest than 
more than 90% of th e central stars in PNe might 
have companions (e.g. iMoe &: De Marcoll2006l ). 

Since the early observation of the PN NGC 2438, which 
is surrounded by stars of the open cluster (OC) M46 
(NGC 2437), photometric and spectroscopic methods have 
been de veloped to test ph ysical associations between such 
objects. IZiz Novsk-^ lll975h published a list of 10 OCs with 
nearby projected PNe. Maiaess. Turner fc Lanel (|2007l ) pro- 
vided a list of 13 possible physical associations and an ad- 
ditional list of 17 angular coincidences {AR < 15'). They 
also present physical criteria for membership based on sim- 
ilarities in velocity, reddening, and the ratio of estimated 
distances. 

aess. Turner &: Land l|2007h reviewed theoretical 
and observational aspects of stellar evolution to impose lim- 
iting parameters for PNe to be OC members. Many co- 
incid ences can be disc arded if the PN is a bulge mem- 
ber (|Ziz Novskvl Il975l ). as indicated by position along 
the disk and radial velocities. OCs younger than 28Myr 



2 C. Bonatto, E. Bica and J.F.C. Santos Jr. 



produce type II SN. The short-hved PN phase varies 
from 10^ to 10^ years, according to the progenitor mass 
ijKoppen fc Ackeji200a ISchonberner fc B15ckeilll99 6'). and 
does not favour associations. PN e observed in globular clus- 
ters fe.g. lKoppen fc Ackeij|2000l ) imply that their turnoffs at 
< 1 Mq are below the lower mass limit for PN progenitors. 
This invokes mass t ransfer in binary systems t o produce PN 
in globular clusters. [Moe fc De Marcoi l|2006l ), ISokeil (|2006h . 
and ,Ziilstra (,2007l 'l presented a scenario where most PNe 
arise from binaries. Consequently, care is necessary because 
the cluster age may not correspond to the PN mass progen- 
itor. 

Amo ng the list of 30 PN / OC coincidences pro- 
vided by iMaiaess. Turner fc Land l|2007l ) , we selected the 
PN/OC pairs NGC 2818/NGC 2818A and NGC2438/M46 
(NGC2437). We also include PK 6-h2.5/NGC 6469 and 
PK 167-0.1, which is projected in the field of an as-yet 
unknown star cluster. We name this low-contrast cluster 
New Cluster 1. In the present paper we derive accurate dis- 
tances to these OCs, which in turn, can be used to investi- 
gate the possible associations. 

Besides the investigation on possible PN/OC associa- 
tions, the present work also serves to provide reliable as- 
trophysical parameters of scarcely studied OCs, as well as 
to characterise a new one. Such data are important both 
to studies of the disc structure and to constrain theories 
of molecular cloud fragmentation, star formation, as well 
as stellar and dynamical evolutions. The present work em- 
ploys near-IR J, H, and Ka photometry obtained from 
the 2MAS^| Point Source Catalogue (PSC). The spatial 
and photometric uniformity of 2MASS, which allow ex- 
traction of large surrounding fields that provide high star- 
count statistics, have been important to derive cluster pa- 
rameters and probe the nature of stellar overdensities (e.g. 
iBica. Bonatto fc CamargollioOTl 'l. 

To this purpose we have developed quantitative tools 
to disentangle cluster and field stars in CMDs, in partic- 
ular two different kinds of filters. Basically we apply (i) 
field-star decontamination to uncover cluster evolutionary 
sequences from the field, which is important to derive red- 
dening, age, and distance from the Sun, and (ii) colour- 
magnitude filters, which proved to be essential for build- 
ing intrinsic stellar radial density profiles (RDPs), as well 
as luminosity and mass functions. In particular, field-star 
decontamination constrains m ore the age and dista nce, es- 
pecially for low- latitude OCs l|Bonatto et al.|[2006ah . These 
tools were applied to OCs and embedded clusters to enhance 
main-sequenc e (MS) and/or p r e-MS sequences with respect 
to the field (iBonatto fc Bica 2006; Bonatto ct al. 2006^; 
Bonatto. Santos Jr. fc Bica 2006; Bica fc Bonatto 2005 



Santos Jr.. Bonatto fc Bical l2005l : iBonatto fc Bical I2005D . tent with = 2660 ± 830 pc (IZhani 



mass function (MF) slopes, have be en derived from colour - 
magnitude-filtered photometry (e.g. IBonatto fc Bicall2005h . 

This paper is organised as follows. Sect. [2] contains ba- 
sic properties and reviews literature data (when available) 
on the possible PN/OC associations. In Sect. |3] we present 
the 2MASS photometry, build CMDs, discuss the field-star 
decontamination, and derive cluster fundamental parame- 
ters. Sect. |4] describes cluster structure by means of stellar 
RDPs and mass density profiles (MDPs). In Sect. O MPs are 
built and cluster masses are estimated. In Sect. [6] aspects re- 
lated to the structure of the present objects are discussed. In 
Sect. [7] we discuss membership considerations. Concluding 
remarks are given in Sect.|8l 



2 THE PRESENT POSSIBLE PN/OC 
ASSOCIATIONS 

In Fig. [1] we show optical XDS^ images of 
NGC 2818/NGC 2818A (left panel, R band) and 
NGC2438/M46 (right panel, DSS B band). Fig. (3 
contains a DSS B image of PK 6-H2.5/NGC 6469 (left panel) 
and an XDSS R image of PK 167-0. 1/New Cluster 1 (right 
panel). 

Table [1] provides fundamental data on the OCs. 
webdaQ coordinates are in cols. 2 and 3; col. 4 gives the 
2MASS extraction radii (Sect.|3|. However, since the RDPs 
(Sect. [4]) built with the 2MASS coordinates presented a dip 
at the centre, new coordinates were searched to maximise 
the central density of stars. The optimised central coordi- 
nates and the corresponding Galactic longitude and latitude 
are given in Cols. 5 to 8 of Table [T] Age, central reddening, 
distance from the Sun and Galactocentric distance based on 
2MASS data (Sect.[3j]) are given in Cols. 9 to 12. Additional 
designations are in col. 13. 

Positional data of the PNe are given in Table [2] where 
we also include the angular separation with respect to the 
OC centre and its relation with the OCs core and limiting 
radii (Sect. [7]). Coordinates are taken from SIMBAeQ. 

Additional information on the possible pairs, related to 
the present work, are summarised below. 

2.1 The pair NGC 2818/NGC 2818A 

Since NGC 2818 refers to the original description of 
the PN, we adopt NGC28 1 8A a s the OC designa- 
tion. IMaiaess. Turner fc Lanj (|2007l) classify t his pa ir as 
a suspected physical association. iPedreroj (|l989l ) de- 
rived a distance from the Sun = 2300 pc and 
a reddening E{B — V) — 0.18, which are c onsis- 

" 1995I ) 



They were useful in the analysis of faint a nd/or distant OCs 
jBica. Bonatto fc D utra 2003; Bica, Bon atto fc Dutrall2004l : 
iBica fc Bonatto! [2OO5I : IBica. Bonatto fc BlumberdbOOd ). In 
addition, more constrained structural parameters, such as 
core and limiting radii {Rc and Ri, respectively), and 



E[B-V) = 0.28 ± 0.15 (|Tvlenda et all ll992D for the 
planetary nebula. Although earlier studies pointed to a 
similar radial ve l ocity between cluster stars and the PN, 
iMermilliod et al.l l|200ll ) give Vr = 20.7 ± 0.3 km s"^ for 
15 red giants of N GC2818A, while for the PN Vr = 
-0.9 ± 2.9kms-i (|Durand. Acker fc Ziilstral 119981 ) and 



1 The Two Micron All Sky Survey, All Sky 

data release fskrutskic ct al. 1997h . available at 
\http://www.ipac.caltech.edu/2mass/releases/allsky/\ 



^ Extracted from the Canadian Astronomy Data Centre 

(CADC), at http:// cadcwww. d ao.nrc. ca/ ' 

obswww.univie.ac.at/webda - \Mermilliod & Paunzerl \200A ) 
* http:/ /simbad. u-starsbg.fr/simbad/ 
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Figure 1. Left panel: 10' X 10' XDSS R image of the OC NGC 2818A and the PN NGC 2818. Rigiit panel: 20' X 20' XDSS R image of 
the OC M 46 and the PN NGC 2438. Images centred on the 2MASS coordinates (cols. 5 and 6 of Table [TJ. North to the top and East to 
the left. 




Figure 2. Left panel: 15' X 15' DSS B image of the OC NGC 6469 and the PN PK 6+2.5. The nearly stellar PN is located close to the 
lower-right corner (see Equatorial coordinates in Table |2j. Right panel: 10' X 10' XDSS R image of the OC New Cluster 1 and the PN 
PK 167-0.1. 



Vr = — 1 =fc 3kms ^ ( Meatheringham. Wood fc Faulkneij 



1 19881 ). iMermilliod et all (|200lh " favour a projection efTeci, 



As for NGC 2818A, a vailable age estimates are 1 Gyr 



l|Mermill iod et al.' '200lJ ) 930 Myr (iLata et all |2002|) and 
794 M yr (Tadross ct al J 120021 '). Additionallv. iTadross et al.1 
l|2002l ) provide for NGC2818A = 2.9 kpc, cluster 

mass m « 288 M0 and the limiting radius Ri — 3.9 pc. 



IMermilliod et al.l (|200lJ ) and iLata et al] ('2002') de- 
rived E{B — V) = 0.18, in agreeme nt with the 
E{B -V) = 0.20 of lTadross et alj(|2002l '). 
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Table 1. General data on the clusters 



Cluster 


a(2000) 


5(2000) 


^cxt 


q(2000) 


5(2000) 


a 


fa 


Ago 


E{B - V) 




dQ 






«GC 




Alternative names 




(hms) 


(°"') 


(') 


(hms) 


C'") 


(°) 


(°) 


(Myr) 






(kpo) 






(kpc) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(S) 


(9) 


(10) 




(11) 






(12) 




(13) 


NGC 2818A 


09:16:01 


-36:37:30 


80 


09:16:07.7 


-36:38:09.6 


262.0 


+ 8.6 


1000 ± 100 


0.10 ± 0.03 


2 


8 + 0. 


1 


8. 


1 + 


.2 


ESO 372- 13, Hen 2-23, PK 261 + 08 


M 46 


07:41:46 


-14:48:36 


80 


07:41:41.1 


-14:51:45.0 


231.9 


+ 4.0 


250 ± 50 


0.10 ± 0.03 


1 


5 + 0. 


1 


8. 


3 + 


.2 


NGC 2437,Mel-75,Cr 159,OC1-601 


NGC 6469 


17:53:13 


-22:19:11 


45 


17:53:03.4 


-22:18:14.4 


6.55 


+ 1.97 


250 ± 50 


0.58 ± 0.06 


1 


1+0. 


1 


6. 


1 + 


.2 


Mel- 182, Cr 353, OCl- 21 ,ESO 589 SC 18 


New Cluster 1 






60 


05:06:20 


+ 39:09:50.0 


167.0 


-1.0 


1000 ± 100 


0.29 ± 0.03 


1 


7 + 0. 


1 


8. 


9 + 


.2 
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Rq = 7.2 kpc 


(tiica et al.||200a). 































2.2 The pair NGC 2438/M46 (NGC 2437) 

Based on the available evidence, aess. Turner fc Land 
l|2007l ) do not exclude the physical association. The distance 
from the Sun and reddening of the PN are = 1775±630 pc 
llZhang|ll995l ) and E{B ~ V) = 0.17 ± 0.08 (|Tvlenda et all 
I I992I ), respectively. The isochrone fit w ith 2MASS photome- 
try bv lMaiaess. Turner fc Land (|2007l ) provided for the OC 
M46 d© = 1700±250pc, E{B -V) = 0.13±0.05, and age = 
22 X lO'^ Myr. Early determinatio ns of radial velocity indi- 
cated differences of « 30kms~^ (|Maiaess. Turner &: Land 
|2007| . and references therein), which suggest spatial coinci- 
dence. However, the most recent s tudy on radial velocity o f 
cluster stars and planetary nebula l|Pauls fc Kohouteklll996l ) 
shows comparable values, reinforcing the possible physical 
association. Additional available data on the OC M 46 are 
age = 250 Myr, dg = l.Skpc core Rc = 2.27 ± O.lSpc 
and l imiting radii -Ri = 11.6 pc (|Nilakshi. Pandev. fc MohanI 
|2002| ): E{B~V) = 0.15, age = 144 Myr, dp = 1.4 kpc, 
i?c = 7.2', and _Ri = 23' (Kharchenko et al.''2005'); i?c = 
4.2±0.3pc and Ri = 11 pc (,Sharma et aI..,2006 V 



2.3 The pair PK 6+2. 5/NGC 6469 

The PN is also known as Ml-31, Ve3-59, ESO 589-PN 16, 
and PNG 006.4-1-02.0. No distance to th e PN is avail- 
able, but iDurand. Acker fc Ziilstral (119981 ) provide Vr = 
68.8 ± 1.8kms"^ The PN is projected on the bulge, 
which decreases the pos sibility of a physical association. 
iKharchenko et al.l (|2005l ) obtained for the OC NGC 6469 
age = 230 Myr, E{B -V) = 0.30, do = 0.55 kpc, R^ = 5.4', 
and _Ri = 9'. 



2.4 The pair PK 167-0. 1/New Cluster 1 

The PN is also known as A 55 7 and PNG 167.0-00.9. 
The available distance and radial velocity estimates are 
da = 1780 pc ([Ph illips' 2004) and Vr = 58.2 ± 6.5 km s'^ 
(|Durand. Acker fc Ziilstra ,1998j ), respectively. The OC 
New Cluster 1 was found by one of us (E.B.) on Digitised 
Sky Survey images, and is analysed for the first time. 



Table 2. PN eoordinates and angular separation with the OCs 



a(2000) 


(5(2000) 


e 


b 


(hms) 


C") 


n 


n 


(2) 


(3) 


(4) 


(5) 



AiJ fc fi 
{') 

(6) (7) (8) 



NGC 2818 09:16:01.7 -36:37:38.8 261.98 +8.58 1.6 0.8 0.1 

NGC 2438 07:41:51.4 -14:43:54.9 231.80 +4.12 8.2 1.8 0.3 

PK6+2.5 17:52:41.4 -22:21:56.8 6.45 +2.01 6.6 6.0 0.5 

PK 167-0.1 05:06:38.4 +39:08:8.6 167.04 -0.97 4.9 10.5 1.4 

Table Notes. Cols. 2-5: coordinates from SIMBAD. Col. 6: angu- 
lar separation of the PN to the optimised cluster centre. Cols. 7 
and 8: ratio of the PN's angular separation with the core {fc = 
AR/Rc) and limiting (/; = AR/Ri) radii (Sect. lU, respectively. 

(cols. 5 and 6 of Table[T} using VizieFilf|. Previous analyses of 
OCs in different environments (Sect, [ij have shown that as 
long as no other populous cluster is present in the field, and 
differential absorption is not prohibitive, wide extraction ar- 
eas provide the required statistics for a consistent field-star 
characterisation in terms of colour and luminosity distribu- 
tion. Thus, the 2MASS extraction radii (col. 4 of Table [TJ 
are significantly larger than the limiting radii (Sect. [4] and 
col. 7 of Table|4]) of the present objects. For decontamination 
purposes, comparison fields were selected within wide rings 
centred on the cluster coordinates and beyond their limiting 
radii. As a photometric quality constraint, 2MASS extrac- 
tions were restricted to stars with magnitudes (i) brighter 
than those of the 99.9% Point Source Catalogue complete- 
ness limil[f|in the cluster direction, and (ii) with errors in J, 
H, and Ks smaller than 0.2 mag. The 99.9% completeness 
limits are different for each cluster, varying with Galactic 
coordinates. The fraction of stars with J, H, and Ks uncer- 
tainties smaller than 0.06 mag is ^ 75% (NGC 6469), ^ 94% 
(New Cluster 1), ^ 73% (NGC2818A), and ^ 93% (M46). A 
typical distribution of uncertainties as a function of 
magnitude, for objects projected to wards the central 
parts o f the Galaxy, can be found in lBonatto Bical 
(|2007bh . Reddening transformations use the relations 
Aj/Av = 0.276, Ah/A v = 0.176, AkJAv = 0-118, and 
Aj = 2.76 X E{J - H) (|Dutra. Santiago fc Bical |2002| ). for 
a constant total-to-selective absorption ratio Rv ~ 3.1. 



3 2MASS PHOTOMETRY 



J, H, and Ks 2MASS photometry was extracted in circu- 
lar fields centred on the optimised coordinates of the objects 



° http://vizier.u-strasbg.fr/viz-bin/VizieR?-source= II/246\ 

^ Following the Level 1 Requirement, according to 

\http://wwwApac.caltech.edu/2mass/releases/allsky/doc/sec6-5al.htrni\ 
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Figure 3. Left: 2MASS CMDs extracted from the R < 5' region of NGC2818A. Top panels: observed photometry with the colours 
J X {J — H) (left) and J X {J — Ks) (right). Middle: equal-area extraction of the comparison field. Besides some contamination of disk 
stars, a populous MS and a giant clump show up. Bottom panels: decontaminated CMDs set with the 1 Gyr Padova isochrone (solid 
line). The colour-magnitude filter used to isolate cluster MS/evolved stars is shown as a shaded region. Right: Same for M46, except for 
the 250 Myr Padova isochrone. 



3.1 Field-star decontaminated CMDs 

J X {J — H) and J x {J — Ks) CMDs of central extrac- 
tions of the clusters are shown in Figs. |3] and |3] Features 
present in the central CMDs and those in the respective 
comparison field (top and middle panels), show that field 
stars contribute in varying proportions to the CMDs, espe- 
cially for the bulge-projected OC NGC6469 (Fig.gll. Never- 
theless, when compared to the equal-area comparison field 
extractions, cluster-like sequences are suggested, especially 
for M 46, NGC2818A, and NGC6469 (the blue sequence). 

To objectively quantify the field-star contamination in 
the CMDs we a p ply th e statistical algorithm described in 
iBonatto fc Bical (|2007bl ). It measures the relative number 
densities of probable field and cluster stars in cubic CMD 
cells whose axes correspond to the magnitude J and the 
colours (J — H) and (J — Ks). These are the 2MASS colours 
that provide the maximum varianc e among CMD sequences 
for O Cs of different ages (e.g. iBonatto. Bica &: Girardj 
|2004| ). The algorithm (i) divides the full range of CMD mag- 
nitude and colours into a 3D grid, (ii) computes the expected 
number density of field stars in each cell based on the num- 
ber of comparison field stars with similar magnitude and 
colours as those in the cell, and (Hi) subtracts the expected 



number of field stars from each cell. By construction, the al- 
gorithm is sensitive to local variati ons of field-star contam i- 
nation with colour and magnitude (|Bonatto fc Bicall2007bl ). 
Typical cell dimensions are AJ = 0.5, and A( J — H) = 
A( J — Ks) — 0.25, which are large enough to allow sufficient 
star-count statistics in individual cells and small enough to 
preserve the morphology of different CMD evolutionary se- 
quences. As comparison field we use wide rings extracted 
from the region Ri < R < 7?oxt around the cluster cen- 
tre to obtain representative background star-count statistics, 
where Ri is the limiting radius (Sect.Q. The comparison 
fields effectively used are located within R — 30' — 80' 
(NGC2818A), R = 40' - 80' (M46), R = 20' - 45' 
(NGC 6469), and R = 10' - 60' (New Cluster 1). In all 
cases, the inner boundary of the comparison field lies 
beyond the probable tidal radius (Sect. |4]). We em- 
phasise that the equal-area extractions shown in the middle 
panels of Figs. [3] and |4] serve only for visual comparison pur- 
poses. Actually, the decontamination process is carried out 
with the large surrounding area as described above. Further 
details on the algorithm, including discussions on subtrac- 
tion effi ciency and limitations, are given in iBonatto fc Bical 
(|2007bl ). 

The decontaminated CMDs are shown in the bot- 
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(J-H) (J-K3) (J-H) (J-K3) 

Figure 4. Same as Fig. |3] for NGC 6469 (left) and New Cluster 1 (right). Isoclirones used are the 250 Myr (NGC 6469) and 1 Gyr 
(New Cluster 1) ones. Despite the heavy contamination by bulge stars, a blue sequence shows up in the CMD of NGC 6469. New Cluster 1 
is a poorly-populated I AC. 



torn panels of Figs. |3] and |4] As expected, most 
of the disk contamination in NGC2818A, M46 and 
New Cluster 1 is removed, leaving stellar sequences typical 
of evolved/intermediate-age OCs. The centrally-projected 
NGC 6469, on the other hand, is so heavily contaminated by 
bulge stars that the algorithm does not subtract all of the 
field stars. The residual bulge component is probably due to 
differential reddening in the field, since the expected Pois- 
son residuals should be smaller. In fact, the residual bulge 
component (bottom panel) contains 398 stars, and the ob- 
served (top panel) contains 1893. Since the cluster sequence 
is significantly bluer than the bulge stars, both sequences 
can be unambiguously separated. We note that in all cases, 
essentially the same CMD features show up in both colours. 
New Cluster 1 has evidence of MS-depletion (Fig. |4]), sug- 
gesting advanced dynamical evolution. 

NGC2818A and M46 have easy-to-detect cluster 
CMD sequences and, to some extent, the same ap- 
plies to NGC 6469. However, the non-populous nature of 
New Cluster 1 requires additional statistical evidence. To 
this end we present in Table |3] the full statistics of the 
decontaminated sequences and field stars, by magnitude 
bins. As a caveat we note that it is more meaningful to 
work with isolated cluster sequences instead of the full 



photometric sample. The bulge residual in the CMD of 
NGC 6469, for instance, contains more stars per magni- 
tude bin than the assumed cluster sequence, which would 
mask the decontamination statistics. In this sense, we 
first isolate the cluster sequences by means of appropriate 
colour-magnitude filters, which are used to exclude stars 
with colours different from those of the assumed cluster 
sequence. They are wide enough to accommodate clus- 
ter MS and evolved star colour distributions, allowing for 
the la photometric uncertainties. Colour-magnitude filter 
widths should also account for formation or dynamical 
evolution-related effects, such as enhanced fractions of bi- 
naries (and other multiple systems) towards the central 
parts of clusters, since such s ystems tend to widen the 



MS ( e.g.lBonatto fc Bica 2007b: B onatto. Bica fc Santos Jr] 
i200a iHurlev fc Tout, 1998. : Kerber et al.ll2002l ). The colour- 
magnitude filters for the present objects are shown in the 
bottom panels of Figs. [3] and |4] 

Statistically relevant parameters to characterise the na- 
ture of a star cluster are: (i) Nio- which, for a given mag- 
nitude bin, corresponds to the ratio of the decontaminated 
number of stars to the Icr Poisson fluctuation of the number 
of observed stars, (ii) crps, which is related to the probability 
that the decontaminated stars result from the normal star 
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count fluctuation in the comparison field and, (in) FSunif, 
which measures the star-count uniformity of the comparison 
field. See below for the precise definition of these parameters. 
Properties of Nio-, crps, and FS unif, measured in PCs and 
field fl uctuations are discussed in lBica. Bonatto fc Camargol 
l|2007h . Table E] also provides integrated values of the above 
parameters, which correspond to the full magnitude range 
spanned by each OC. The spatial regions considered here 
are those sampled by the CMDs shown in the top panels of 
Figs. E] and H 

CMDs of star clusters should have integr ated Nio- val- 
ues s ignificantly larger than 1 (Bica, Bonatt o &: Camargol 
l2007l 'l. Indeed, this condition is met by NGC2818A, M46, 
and NGC6469. New Cluster 1, on the other hand, occurs at 
the ~ 2.7(7 level, close to the lower limit obs erved in central 
directions (jBica. Bonatto fc Camargol I2OO7I ). Values higher 
than 3(T occur as well for the Nio- measured in magnitude 
bins. Because of the small number of bright stars, we point 
out that this analysis should be basically considered for 
J > 10. As a further test of the statistical significance of the 
above results we investigate star count properties of the field 
stars. First, the comparison field is divided in 8 sectors, each 
with 45° of opening angle around the cluster centre. Next, 
we compute the parameter ctfs , which is the 1 a Poisson 
fiuctuation around the mean of the star counts measured 
in the 8 sectors of the comparison field (corrected for the 
different areas of the sectors and cluster extraction). For a 
spatially uniform comparison field, aps should be very small. 
In this context, star clusters should have the probable num- 
ber of member stars (Nci) higher than ~ 3(tfs, to minimise 
the probability that Nd results from the normal fiuctuation 
of a non-uniform comparison field. Again, this condition is 
fully satisfied, in some cases reaching the level Nd ~ 10 ctfs. 
The ratio decreases somewhat for New Cluster 1, but still is 
higher than ~ 3, probably because it is almost projected 
against the anti-centre. Finally, we also provide in Table [3] 
the parameter FSunif. For a given magnitude bin we first 
compute the average number of stars over all sectors (TV) and 
the corresponding la fiuctuation CT(]v>; thus, FSunif is de- 
fined as FSunif ~ CT(N)/(N). Non uniformities such as heavy 
differential reddening should result in high values of FSunif- 

The three statistical tests applied to the present sample, 
i.e. (i) the decontamination algorithm, (ii) the integrated 
and per magnitude Nio- parameter, and (in) the ratio of Nd 
to Ctfs, produce consistent results. As for New Cluster 1, the 
number of decontaminated stars in each magnitude bin is 
significantly larger than what could be expected from field- 
star fiuctuations. Besides, its field is rather uniform (Ta- 
ble[3]). We conclude that New Cluster 1 is a poorly-populated 
Gyr-class OC. 



l2004h . The decontaminated CMD morphologies (bottom 
panels of Figs. [3] and (4)1 provide enough constraints to de- 
rive reliable cluster fundamental parameters, which are given 
in Table [1] The isochrone fits to the decontaminated 
cluster CMDs are shown in the bottom panels of 
Figs. [3] and [H 

From the isochrone fit to the CMD of 
NGC 2818 A (Fig. [a]) we derive the age 1.0 ±0.1 Gyr, red- 
dening E{J - H) = 0.03±0.01 that converts to E{B - V) = 
0.10 ± 0.02 and Ay = 0.31 ± 0. 07, and dp, = 2.8 d= 0.1 kpc. 
We adopt Rq = 7.2 ± 0.3 kpc (i Bica et ahlbOOeh as the 
Sun's distance to the Galactic centre to compute 
the OC's Galactocentric distances. The latter value 
was derived by means of the globular cluster spatial 
distribution. Besides, other recent studies gave sim- 
ilar results, e.g. Rq — 7.2 ±0.9 kpc (Eisenhauer et aU 
[2003), Rq = 7.62 ± 0.32 kpc (Eisenhauer et al. 200^ 
and Rq = 7.52 ±0.10 kpc (jNishivama et al.l 12006). 
with different approaches. Thus, for Rq = 7.2 kpc, 
the Galactocentric distance of NGC2818A is Rqc = 8.1 ± 
0.2 kpc. NGC 2818A lies f» 0.9 kpc outside the Solar 
circle. Within the uncertainties, t he present age and red- 
dening values ag r ee wi th those of iLata et al] (|2002l ) and 
iMermilliod et all (12001 ). The di s tance from the Sun is the 
same as that in iTadross et all (|2002l ). The mass at the 
turnoff (TO) of NGC2818A is mxo = 2.1 M©. 

Parameters of M 46 are the age 250 ± 50 Myr, 
E[J~H) = 0.03 ± 0.01, E{B~V) = 0.10 ± 0.02, 
and Av = 0.31 ± 0.07, cLq = 1.5 ± 0.1 kpc, and 
^GC = 8.3 ± 0.2 kpc, ~ 1.1 kpc outside the So- 
lar circle. Within the uncertainties, these parameters 



agree with those of Maiaess. Turner fc Land (i20P7) and 
iNilakshi. Pandev. fc MohanI ll20o3r th e pres ent age is about 
twice as that of iKharchenko et al.l |200^. TO mass is 
mTO = 3.5 Mq. 

For the bulge-projected NGC 6469 we derive the age 
250 ± 50 Myr, E(J - _ff) = 0.18 ± 0.02, E{B - V) = 0.58 ± 
0.05 and Av = 1.8±0.2, = 1.1±0.1 kpc, and 7?gc = 6.1± 
0.2 kpc, ~ 1.1 kpc inside t he Solar circle. T he present 
age agrees with that of Kha rchenko et al] (|2005l ). but they 
found about twice the reddening and half the distance from 
the Sun. TO mass is mxo = 3.5 Mq. 

Finally, parameters of New Cluster 1, derived for the 
first time, are the age l.OiO.l Gyr, E{J - H) = 0.09±0.01, 
E{B-V) = 0.29 ± 0.03, and Av = 0.9 ± 0.1, = 
1.7 ±0.1 kpc, and Rgc = 8.9 ± 0.2 kpc, « 1.7 kpc outside 
the Solar circle. TO mass is mxo = 2.1 Mq. 

Except for the bulge-projected NGC 6469, which is lo- 
cated ~ 1 kpc inside the Solar circle, the remaining OCs are 
more than 1 kpc outside it. 



3.2 Cluster age, reddening, and distance 

Fundamental parameters for the present clusters are derived 
from CMD fits wi th solar-metallicity Padova isochrones 
(jCirardi et all 120021 ) computed with the 2MASS J, H, and 
Ks filterfl These is ochrones are very sim ilar to the Johnson- 
Kron-Cousins (e.g. iBessel fc Brett 19881) ones, with differ- 
ences of at most 0.01 in (J — H) ( Bonatto. Bica fc Girardl 



''' I http://stev. oapd. inaf. it/r^lgirardi/ cgi- bin/ cmd\ 



4 CLUSTER STRUCTURE 

Structural parameters are derived by means of RDPs, de- 
fined as the projected radial distribution of the number den- 
sity of stars around the cluster centre. 

Star clusters usually have RDPs that follow some well- 
defined analytical profile. The most often used ar e the single 
mass, modified isotherma l sphere of Kind (|l966al ). the mod- 
ified isotherm al sphere oflWilsonI (Il975|), an d the power law 
with a core of lElson. Fall fc Freeman " (|1987l ). Each function 
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Table 3. Statistics of the field-star decontamination in magnitude bins 



Colour-magnitude filtered photometry 



NGC2818A M46 NGC6469 New Cluster 1 



AJ 


cps 


Nobs 


Nci 


Ni,, 


FSunif 


0"FS 


Nobs 


Nci 


Ni„ 


FSunif 




Nobs 


Nci 


Ni,, 


FSunif 


0"FS 


Nobs 


Nci 


Ni,, 


FSunif 


8-9 












0.11 


3 


3 


1.7 


0.28 


0.14 


1 


1 


1.0 


0.37 












9-10 


0.14 


1 


1 


1.0 


0.22 


0.25 


2 


2 


1.4 


0.39 


0.31 


3 


3 


1.7 


0.45 


0.29 


2 


2 


1.4 


0.45 


10-11 


0.34 


2 


2 


1.4 


0.10 


0.30 


12 


12 


3.5 


0.53 


0.74 


5 


5 


2,2 


0.46 












11-12 


0.34 


18 


17 


4.0 


0.07 


0.31 


22 


22 


4.5 


0.26 


1.92 


14 


12 


3.2 


0.40 


0.77 


11 


11 


3.3 


0.17 


12-13 


1.36 


20 


12 


2.7 


0.14 


0.28 


43 


43 


6.3 


0.06 


2.83 


36 


25 


4.2 


0.16 


1.07 


7 


5 


1.9 


0.15 


13-14 


1.43 


31 


27 


4.8 


0.25 


1.16 


72 


72 


6.4 


0.05 


7.01 


67 


27 


3.3 


0.14 


2.09 


19 


3 


0.7 


0.11 


14-15 


3.47 


78 


63 


7.1 


0.18 


3.18 


85 


85 


5.4 


0.08 












3.35 


40 


12 


1.9 


0.08 


15-16 


4.51 


81 


36 


4.0 


0.09 
































All 


10.95 


231 


158 


9.9 


0.12 


4.97 


238 


182 


11.6 


0.07 


9.20 


126 


73 


5.9 


0.12 


6.42 


79 


33 


2.7 


0.09 



Table Notes. This table provides, for each magnitude bin (AJ), the 1 a Poisson fluctuation ((Jfs) around the mean, with respect to the 
star counts measured in the 8 sectors of the comparison field, the number of observed stars (Nobs) within the spatial region sampled in the 
CMDs shown in the top panels of Figs. [3] and |4] the respective number of probable member stars (Nd) according to the decontamination 
algorithm, the Ni^ parameter, and the field-star uniformity parameter. The statistical significance of Nd is refiected in its ratio with the 
1(7 Poisson fluctuation of f^ohs (Nict) and with crpg. The bottom line corresponds to the full magnitude range. 



is characterised by different parameters that are somehow 
related to cluster structure. However, because the error bars 
in the present RDPs are significant (see below), we decided 
to use the analytical profile a (J?) = atg + <jo/{1 + (R/ Rc)^), 
where cri,g is the residual background density, ao is the cen- 
tral density of stars, and Rc is t he cor e radiu s. This function 
is similar to that introduced bv lKind (|l962l l to describe the 
surface brightness profiles in the central parts of globular 
clusters. 

In all cases we build the stellar RDPs with colour- 
magnitude filtered photometry (Sect. [ST). However, residual 
field stars with colours similar to those of the cluster are ex- 
pected to remain inside the colour-magnitude filter region. 
They affect the intrinsic stellar RDF in a degree that de- 
pends on the relative densities of field and cluster stars. The 
contribution of the residual contamination to the observed 
RDF is statistically taken into account by means of the com- 
parison field. In practical terms, the use of colour-magnitude 
filters in cluster sequences enhances the contrast of the RDF 
with respect to the background level, espe cially for objects 
in dense fields (e.g. iBonatto fc Bicair2007bl ). 

To avoid oversampling near the centre and undersam- 
pling at large radii, RDFs are built by counting stars in 
rings of increasing width with distance from the cen- 
tre. A typical distribution of ring widths would be 
Ai? = 0.5, 1, 2, 5, and 10', respectively for 0' R < 1', 
1' s; i? < 4', 4' ^ < 10', 10' ^ 7? < 30', and R > 30'. 
However, the number and width of the rings are adjusted 
to produce RDFs with adequate spatial resolution and as 
small as possible la Foisson errors. The residual background 
level of each RDF corresponds to the average number of 
colour-magnitude filtered stars measured in the comparison 
field. The R coordinate (and respective uncertainty) of each 
ring corresponds to the average position and standard devi- 
ation of the stars inside the ring. 

The resulting RDFs of the present star clusters are 
given in Fig. O For absolute comparison between clusters 
the radius scale was converted to parsecs and the num- 
ber density of stars to stars pc~^ for the distances derived 



in Sect. 13.21 Besides the RDFs resulting from the colour- 
magnitude filters, we also show, for illustrative purposes, 
those produced with the observed (raw) photometry. In all 
cases, minimisation of the number of non-cluster stars by 
the colour-magnitude filter resulted in RDFs with a signifi- 
cantly higher contrast with respect to the background. Fits 
of the King-like profile were performed with a non-linear 
least-squares fit routine that uses errors as weights. To min- 
imise degrees of freedom, ao and Rc were derived from the 
RDF fit, while atg is measured in the respective compari- 
son field. These values are given in Table 3] and the best- 
fit solutions are superimposed on the colour-magnitude fil- 
tered RDFs (Fig. 0. Because of the 2MASS photometric 
limit, which in most cases corresponds to a cutoff for stars 
brighter than J « 15, ao should be taken as a lower limit 
to the actual central number density. The adopted King-like 
function describes well the RDFs throughout the full radii 
range, within uncertainties. 

We also estimate the cluster limiting radius and un- 
certainty by visually comparing the RDF level (taking into 
account fluctuations) with the background. In this sense, Ri 
corresponds to the distance from the cluster centre where 
RDF and background become statistically in distinguishable 
from each other (e.g. iBonatto fc Bicall2005l . and references 
therein). For practical purposes, most of the cluster stars can 
be considered to be contained within Ri. Note that Ri should 
not be mistaken for the tidal radius. For instance, in popu- 
lous and relatively high Galactic latitude OCs such as M 26, 
M67, NGC188, and NGC 2477, limiting radii arc a factor 
~ 0. 5 — 0.7 of the respective tidal radii (Bonatto & Bical 
l2005h . The limiting radii of the present objects are given in 
col. 7 of Table |4l Tidal radii are derived from fits of King 
profile to RDFs, which depend on wide surrounding fields 
and adequate Foisson errors. If limiting and tidal radii 
of the present clusters are similarly related as for 
the bright ones, we note that, in all cases, the lower- 
limit of the radial range adopted as comparison field 
(Sect. [sTTj) is located beyond the respective probable 
tidal radius. This, in turn, minimises the probability 
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Table 4. Structural parameters of evolved + MS stars (CM-filtcred photometry) 



RDP MDP 



Cluster 


1' 






Sc 


Rc 


Ri 


CO 


Rc 




(pc) 


(stars pc~^) 


(stars pc~^) 




(pc) 


(pc) 


(Mqpc-2) 


(pc) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 2818A 


0.805 


1.8 ±0.1 


8.8 ±2.2 


5.9 ± 1.2 


1.5 ±0.3 


11 ± 1 


13.4 ± 1.8 


1.5 ±0.2 


M46 


0.442 


4.4 ±0.1 


18.8 ±4.6 


5.1 ± 1.0 


2.0 ±0.4 


12.4 ±0.8 


30.8 ±6.9 


2.0 ±0.3 


NGC 6469 


0.318 


9.8 ±0.2 


35.3 ± 13.0 


4.6 ± 0.5 


0.35 ±0.10 


4.0 ±0.5 


49 ±8 


0.45 ±0.10 


New Cluster 1 


0.493 


3.7 ±0.1 


20.8 ±6.6 


6.6 ± 1.8 


0.23 ±0.06 


1.7 ±0.2 


27.9 ± 1.6 


0.25 ±0.02 



Table Notes. Col. 2: arcmin to parsec scale. King profile is expressed as cr(R) = (Ji,g + ito/(1 + (R/ Rcore)^)- To minimise degrees of 
freedom in RDP fits, crf,g was kept fixed (measured in the respective comparison fields) while ag and Rc were allowed to vary. MDPs are 
background subtracted profiles. Col. 5: cluster/background density contrast {5c = 1 + o-o/o-i,g), measured in CM-filtered RDPs. 



of cluster members at large radii, e.g. in the halo, to 
be considered as field stars by the decontamination 
algorithm. 

Table |4] (col. 5) also provides the density contrast pa- 
rameter Sc = 1 + cTo/o'bg- Since Sc is measured in colour- 
magnitude-filtered RDPs, it does not necessarily correspond 
to the visual contrast produced by observed stellar distri- 
butions in images (Figs. [1] and [2} . NGC 6469, for instance, 
presents a very low contrast in the DSS B image (Fig. [2]) but, 
because most of the non-cluster stars have been excluded by 
the colour-magnitude filter, the corresponding RDP presents 
a relatively high density contrast, Sc « 4.6. The same applies 
to the non-populous OC Nevy Cluster 1. 

Probably because of different methods and data sets, 
the present valu es of Rc and Ri a r e diff erent from those 
in common with iKharchenko et ahl (|2005l ) . The difference, 
especially in Rc, may be attributed to their brighter limits 
l|Kharchenko et al.l20oi ) producing shallower profiles for the 
bulge-contaminated cluster NGC 6469. 



4.1 Mass Density Profiles 

To complete the structural description of the objects we take 
the mass-luminosity (ML) relation derived from isochrone 
fits (Sect. [T2)) to build statistical mass-density profiles. We 
follow the same systematics as that used to build RDPs. 
Instead of computing the number density of stars in rings, 
we now assign each star a mass according to the respective 
ML relation. MDPs are produced by subtracting from the 
mass density in each ring that measured in the comparison 
field. They are shown in Fig. |6l together with the respective 
King-like fits. Likewise RDPs, MDPs are well described by 
King-like profiles. Core radii derived from MDPs (col. 9 of 
Table|4ll agree, at la, with RDP ones (col. 6). 
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Figure 5. Stellar RDPs (filled circles) of the present star clusters 
built with colour-magnitude filtered photometry. Solid line: best- 
fit King-like profile. Horizontal shaded region: offset field stellar 
background level. Gray regions: la King fit uncertainty. The core 
and limiting radii are indicated. Insets: RDPs built with the ob- 
served photometry. Absolute scale is used. 



5 MASS FUNCTIONS AND CLUSTER MASS 

The methods presented in iBonatto fc Bical l|2005l ) (and 
references therein) are used to build the mass func- 
tions (MFs), {(j){m) = ^). We build them with colour- 
magnitude filtered photometry, the three 2MASS bands 
separately, and the ML relations obtained from the re- 
spective Padova isochrones and distances from the Sun 



(Sect. | 3.2p. Further details o n MF construction are given in 
iBica. Bonatto fc Blumbergj (|2006l V The effective magnitude 
range over which MFs are computed is that where clusters 
present an excess of stars with respect to the comparison 
field. In all cases it begins right below the TO and ends at 
a faint magnitude limit brighter than that stipulated by the 
2MASS completeness limit (Sect.EJ. The effective MS stellar 
mass ranges are 1.1 m(MQ) sC 2.1 (NGC2818A), 1.1 ^ 
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Figure 6. Same as Fig. [5]for the mass density profiles. The back- 
ground contribution has been subtracted. 
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Figure 7. Overall mass functions fitted with (f>{ni) oc m~(^+^'. 
Within the uncertainties, the MP slopes are similar to Salpeter's 
IMF. 



m(M0) < 3.5 (M46), 1.3 < m{MQ) < 3.5 (NGC6469), and 
1.1 m{MQ) ^ 2.1 (NewCluster 1). However, we note that 
because of the non-populous nature of NewCluster 1, the 
MF error bars resulted exceedingly large. 

The resulting MFs are shown in Fig. [T] where fits with 



the function 



(ml cx m 



-(i+x) 



are included; MF slopes are 



given in col. 4 of Table [S] The populous nature of M46 
allowed computation of the core and halo MF parameters; 
for NGC2818A and NGC6469, only the overall MF was 
considered. Within uncertainties, th e MFs o f the present 
clusters have slopes similar to that of lSalpete 3 11195511 Initial 
Mass Function (IMF) (x = 1.35). The core MF of M46 
appears flatter than the halo's. 

Table [5] gives parameters of the target clusters mea- 
sured in the CMDs and derived from MFs. The number 
of evolved stars (col. 2) was obtained by integration of the 
(background-subtracted) colour-magnitude filtered luminos- 
ity function for stars brighter than the TO. This number 
multiplied by the stellar mass at the TO yields an estimate 
of the mass stored in evolved stars (col. 3). The number and 
mass of the observed MS stars (cols. 5 and 6, respectively) 
were derived by integrati the MFs over the effective MS mass 
range. 

To estimate the total stellar mass we extrapolate the 
observed MFs down to the H-burning mass limit (0.08 A/q). 
For masses below the presen t detection thr eshold (Table [5]) 
we base the extrapolation on iKroupal (|200ll ) universal IMF, 
in which x = 0.3 ± 0.5 for the range 0.08 < m{MQ) < 0.5 
and X = 1.3±0.3 for 0.5 < m(MQ) < 1.0. When the present 
MF slopes are flatter than or similar (within uncertainties) 



to Kroupa's, we adopt the measured values of x- The to- 
tal (extrapolated MS -|- evolved) values of number, mass, 
surface, and volume densities are given in cols. 7 to 10 of 
Table [i 

The mass stored in the observed MS and evolved stars 
of M 46 amounts to Mobs ~ 1 500 Mq , compatible with the 
populous CMD (Fig. [3]). The present observed mass esti- 
mate for NGC2818A, Mpbs ~ 240 M© agrees with the es- 
timate of iTadross et al.1 (|2002l ). NGC6469, and especially 
NewCluster 1, are low-mass OCs, respectively with Mobs ~ 
120 AIq and Mobs ~ 30 Mq . The extrapolated masses are a 
factor ~ 5 times larger than the observed ones. As a caveat 
we note that the total mass estimates should be taken as 
upper limits, since because of dynamical evolution, signifi- 
cant fractions of the low-mass content may have been lost 
to the field. 



6 COMPARISON WITH NEARBY OCS 

At this point it is interesting to compare the struc- 
tural parameters derived for the present OCs with 
those of a reference sample of nearby OCs with 
ages in the ra nge 70 — 7 000 Myr and masses within 
400 - 5 3OOM0 (|Bonatto fc Bic a 2005). To the original ref- 
erence sample were added t he yo ung OCs NGC6611 
( Bonatto. Santos Jr. fc B"ical |2006| ) and NGC4755 
( Bonatto et al] l2006bl V Clusters are differentiated ac- 
cording to total mass (smaller or larger than IOOOMq). 
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Table 5. Parameters related to cluster mass and stellar population 



Evolved Observed MS Evolved + Extrapolated MS 

Region N* m x ^* "nioba ^* ma p 

(stars) (Mq) (stars) (Mq) (lO^stars) (lO^Afo) (Mqpc-^) (Mq pc-^*) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 

NGC 2818A — MS: 1.1 ^ m{MQ) ^ 2.1 — Age = 1.0 ± 0.1 Gyr 

Core 6±114±3 — 24 ±5 40 ±5 — — — — 

Halo 18 ± 8 38 ± 18 — 90 ± 20 149 ± 20 — — — — 

Overall 25 ± 9 52 ± 18 1.44 dz 0.55 115 ± 20 189 dz 21 41 ± 36 15 ± 8 3.9 ± 2.0 0.26 ±0.14 

M46 — MS: 1.1 ^ m{MQ) ^ 3.5 — Age = 250 ± 50Myr 

Core 9± 2 29 ± 6 1.31 ±0.50 133 ± 15 214 ± 14 17 ± 13 7±3 55 ± 21 21 ± 8 

Halo 22 ±6 68 ± 19 1.60 ±0.22 727 ± 70 1147 ± 56 112 ± 83 43 ± 16 9.0 ± 3.3 0.54 ±0.20 

Overall 32 ± 6 98 ± 20 1.63 ±0.15 860 ± 71 1362 ± 58 135 ± 99 51 ± 18 10.6 ± 3.8 0.64 ± 0.23 

NGC 6469 — MS: 1.3 ^ m(MQ) ^ 3.5 — Age = 250 ± 50Myr 

Core — — — 9±4 15±3 — — — — 

Halo — — — 65 ±31 107 ±28 — — — — 

Overall 1±1 2±2 1.38 ±0.63 74 ± 31 122 ± 29 16 ± 13 6.3 ± 2.8 13 ± 6 2.4 ± 1.0 

New Cluster 1 — MS: 1.1 ^ m(MQ) < 2.1 — Age = 1.0 ±0.1 Gyr 

Core — — — — — — — — — 

Halo — — — — — — — — — 

Overall 7±2 14±4 — 13 ±8 18 ±7 — — — — 



Table Notes. Col. 6: stellar mass stored in the observed MS. Col. 8: mass of the evolved stars added to the MS mass extrapolated to 
0.08 Mq. 



Details on parameter correlation in the reference sample 
are given m iBonatto fc Bical (|2005l ). 

In panels (a) and (b) of Fig. [8] we compare limiting 
and core radii of the present OCs with those of the refer- 
ence sample in terms of cluster age. With respect to both 
kinds of radii, the bulge-projected NGC 6469 and the poorly- 
populated New Cluster 1 appear to be significantly smaller 
than the nearby OCs of similar ages, especially in limit- 
ing radius. Core and limiting radii in the reference sam- 
ple are related by Ri — (8.9 ± 0.3) x iiciro*"'^'' (panel 
(c)), which suggests a similar scaling for both kinds of 
radii, at least for the radii ranges 0.5 < Rc{pc) < 1.5 and 
5 < Ri{pc) < 15. Within uncertainties, NGC2818A, M46, 
NGC 6469, and New Cluster 1 also follow that relation. Fi- 
nally, except for New Cluster 1, the remaining OCs appear 
to follow the trend of increasing limiting radii with Galac- 
tocentric distance (panel (d)). A similar dependence with 
Roc holds as well for Rc, because of the relation implied by 
panel (c). 



7 DISCUSSION 

Besides the spatial coincidence, for a physical association 
between a PN and a star cluster to be considered as highly 
probable, at least the radial velocities, reddening values, and 
distances should be compatible. In the present work we de- 
rive reliable values of reddening, distance from the Sun, core 
and limiting radii for the OCs. The angular separation of the 
PNe NGC 2818, NGC 2438, and PK 6+2.5 with respect to 
the OC centres (Table [2]), corresponds to the halo; PK 167- 
0.1, on the other hand, appears to lie near the border of 



tr' 




0.1 1 6 7 8 9 



Rco,» (Pc) d3c(kpc) 

Figure 8. Relations involving structural parameters of OCs. Cir- 
cles: nearby OCs, including two young ones. Filled symbols: OCs 
more massive than 1 000 Mq . 
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New Cluster 1. In all cases, the mass at the TO of the OC is 
consistent with the presence of a PN. 

Based on the above criteria, the best candidate to 
a physical pair is the PN NGC 2438 with the OC M46. 
They have comparable radial velocities (Sect. 12. ip . and 
the reddening and distance from the Sun for the PN are 
E{B-V) = 0.17 ± 0.08 and d© ~ l.Skpc (Sect. 
while for the OC we derive E(B - V) = 0.10 ± 0.03 and 
= 1.5 ± 0.1 kpc. Besides, the presence of a PN in M46 is 
compatible with the TO mass, rriTO ~ 3.5 Mq. 

Within uncertainties, the present reddening 
E{B ~V)= 0.10 ± 0.03 for the OC NGC 2818A is compat- 
ible with the E{B - V) = 0.18 and E{B - V) = 0.28 ± 0.15 
(Sect, [m for the PN NGC 2818. Distances from the Sun 
are in better agreement, = 2.8 ± 0.1 kpc for NGC 2818A 
and d© = 2.3 - 2.7kpc for NGC 2818 (Sect.[2l]). However, 
the radial velocities do not agree (Sect. I2.1|l . which could 
indicate a projection effect. 

PK 6+2.5 does not appear to be physically associated 
with NGC 6469. If this nearby cluster (Table[l| has a circular 
orbit about the Galactic centre, a small value for its radial 
velocity is expected from its Galactic location. However, the 
relatively high velocity Vr = 68.8 ± 1.8 km s"^ (Sect. [23]) of 
PK 6+2.5 suggests rather a relation to the bulge. 

The distance of 1780 pc of PK 167-0.1 (Sect, i^]) is 
in excellent agreement with that derived for New Cluster 1, 
d© = 1.7±0.1 kpc (Table[l|. This piece of evidence suggests 
physical association, but the cluster radial velocity should 
be determine d for a more robust co mparison. 

Majacss, Turner fc Lan3 (|2007l ) comment that most po- 
tential Planetary Nebulae in OCs are projected against the 
cluster halo. At least in part, this observation can be a 
direct consequence of the relatively high cluster halo/core 
cross section ratio. Ah/ Ac = {Ri — Rc)^/Rc^. Indeed, from 
Fig. |S] (panel c) it follows that the limiting and core radii 
are related by _Ri ~ 9 Rc, which implies that, for a typi- 
cal OC, Ah /Ac ~ 60. In this sense, random distributions of 
(spatially unrelated) PNe and OCs would result in a signif- 
icantly larger number of PNe projected against halos than 
central regions. Fo r the 30 spatial coinciden c es (w ith sepa- 
ration AR < 15''). iMaiaess. Turner fc Land l|2007l ') present 
an estimate of a nuclear radius (r„) and measure the angu- 
lar separation of the PN with respect to the cluster centre. 
Although their r„ is not derived from radial profiles, it prob- 
ably represents the core radius to some degree. For the two 
additional pairs dealt with in the present paper, the ratio of 
the number of cluster halo/core PNe is Nh/Nc ~ 26/6 — 4.3, 
less than 10% of that expected of random, unrelated distri- 
butions. Interestingly, for a King-like cluster characterised 
by the core and limiting radii _Rc and _Ri, the number of 
member stars in the halo is related to that in the core by 
Nh/Nc = In [l + {Ri/Rcf] /ln(2) - 1. Thus, for Ri « 9i?c, 
we obtain Nh/Nc « 5.4, close to the cor responding observed 
ratio for the PNe of iMaiaess. Turner fc Lanel (120071 '). 

We further examine this issue from a dif- 
ferent perspective, sim i lar to that employed by 
iLundstrom &: Stenholml (1984) to study spatial co- 
incidences of WR stars with OCs. We compute the 
ratio of the PN/OC separation to the cluster ap- 
parent radius x ~ A R/Rpc, based on the angu- 
lar diameters given by iDias et al.1 (|2002l '). Accord- 
ing to this definition, x represents the PN sep- 



aration as a function of the apparent cluster ra- 
dius. Next we build the surface density distribu- 
tion of the measured x values, i.e, the number of 
spatial coincidences characterised by x P^r clus- 
ter area, which is shown in Fig. [9l To investigate 
the significance of this distribution we first test 
whether it could result from a random spatial dis- 
tribution of PNe. The 30 PN/OC angular separa- 
tions in iMai aess. Turner &c Land (HoOT.), and the 
2 additional ones (Table [2]) not in common with 
them, are restricted to AR < 55'. For each of the 
32 OCs corresponding to the above spatial coinci- 
dences, we randomly take values of AR uniformly 
distributed in the range 0' < A i? < 55' and com- 
pute X foi" the actual apparent OC radius. The re- 
spective surface density distribution was built for 
lO'' runs, which was subsequently normalised to 32 
objects. The measured curve presents a conspicu- 
ous excess over the simulated one for X 3, which 
might indicate a physical relation. Both curves agree 
for X ^ 3, within the uncertainties. The second test 
follows ILundstrom &; Stenholm (1984): the sign of 
the Galactic latitude (6) of each PN is reversed, but 
keeping the original £. Based on the new PN coor- 
dinates, we then searched for an OC that is closer 
than AR< 55', a condition that was matched for 20 
cases. Within the uncertainties, the corresponding 
surface density distribution, normalised to 32 ob- 
jects (Fig. [9|), coincides with the simulated (random 
spatial distribution of A i? values) one along the full 
range of x- 



The arguments above are based on the appar- 
ent cluster radii of tolas et al.l (|2002l ). Because they 
are basically visually determined, such radii, how- 
ever, are significantly underestimated with respect 
to the limiting radii and, especially, the tidal radii. 
For in st ance, the OCs inlB ica. Bonatto & Blumberd 
(2006'), 'Bon atto &: BicaT (2007b), Bonat to Bical 
(2006), and B onatto et al. (2006b). which are similar 
to those in Maiae ss. Turner &: Land (12007) and were 
studied with the same methods as in the present 
paper, have a limiting radius ~ 2 — 3 times larger 
than the apparent radius given in .Dias et al.l (|2002h . 
Thus, if a similar relation holds for the apparent and 
limiting radii in the present OC sample, the excess 
in the observed surface density (Fig. would cor- 
respond to PN separations AR < Ri « Rt, which 
again suggests a physical relation. 



Despite the many uncertainties associated with the 
above arguments, the small number of targets and the 
uncertain PN distances in particular, it might, nev- 
ertheless, be speculated whether the relatively large frac- 
tion of core PNe results from a physical relation. A cau- 
tionary remark is that the rema ining spatial coincidences in 
[Mai aess. Turner fc Land l|2007l ) should be analysed in detail 
for a better definition of the core and limiting radii, as well 
as the fundamental parameters. 
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O o PN/OCs of the present paper 

■ — • Random distribution of ARp^^,^ 
- ^^ Simulation: reverse b 




Figure 9. Surface density distribution (number of spatial coin- 
cidences characterised by x per cluster area). Empty circles: the 
32 spatial coincidences in this paper. Filled circles: distribution 
simulated with 10® values of A R, randomly taken from the range 
0' — 55', divided by the actual apparent cluster radii. Empty tri- 
angles: distribution of the 20 spatial coincidences that result after 
reversion of the sign of PN's b coordinate. The latter two distrib- 
utons were noma lised to 32 objects. The apparent cluster radii 
l|Dias et al.l l2002h are ~ 2 — 3 times smaller than the limiting 
radii. 



sistent with physical association. In the cases of physical 
association, the turnofT masses are compatible with the oc- 
currence of PNe. It would be important to determine 
radial velocities for members of the new cluster, to 
compare them with that of PK 167-0.1, in order to 
further test the physical association possibility. Be- 
sides, to establish the non-association of PK 6+2.5 
with NGC 6469, radial veocities of member stars are 
necessary. 

From the PN/OC angu lar positions provided by 
iMaiaess. Turner fc Lanj (|2007l ). we estimate that the frac- 
tion of PNe projected against the central region of clus- 
ters, with respect to those in the halo, is significantly higher 
than that expected from spatially unrelated distributions of 
PNe and OCs. Besides, this fraction agrees with that ex- 
pected of clusters characterised by a King-like radial stel- 
lar distribution. In addition, we analysed the distribu- 
tion of the number of spatial coincidences per clus- 
ter area built with the 32 cases dealt with in this 
paper. This curve was compared to that expected 
from a spatially-unrelated distribution of PN/OC 
separations and that resulting from inversion of the 
sign of the b PN coordinates. The observed distri- 
bution presents a conspicuous excess, for PN/OC 
separations smaller than the cluster limiting radii, 
over both simulations, which might indicate phys- 
ical relation. However, we note that these results 
do not take into account the rather uncertain PN 
distances, and radial velocities. Besides, they are 
based on a relatively small number of spatial co- 
incidences. In this context, it is essential to analyse 
in more detail (i.e. with accurate fundamental and 
structural parameters) the wh ole PN/OC sample of 
iMaiaess. Turner &c Land (|2007i ) in order to inves- 
tigate whether the above excess of relatively close 
spatial coincidences can be interpreted as physical 
association. 



8 SUMMARY AND CONCLUSIONS 

In the present study an open cluster was discovered 
(New Cluster 1), which is located at a(J2000) = 05''06™20" 
and 5(J2000) = +39°09'50". Projected on its field, with 
an angular separation of 4.9' with respect to the 
custer centre, is the PN PK 167-0.1. Field-star decon- 
taminated 2MASS photometry was employed in the analysis 
of 4 pairs of PN/OC candidates to physical association. We 
derived accurate ages, distances from the Sun, reddening 
values, and the mass at the turnoff for the clusters. Cluster 
core and limiting radii were determined by means of radial 
density profiles fitted with a King-like function. 

The values of reddening and distance from 
the Sun derived in this work suggest that 
the pairs PN/OC NGC 2818/NGC 2818A and 
NGC2438/M46 are compatible with physical asso- 
ciations. However, the available radial velocity data 
for NGC 2818/NGC 2818 A favour a projection ef- 
fect. The pair PK 6-1-2. 5/NGC 6469 appears to be physically 
unrelated, since the radial velocity of PK 6-1-2.5 suggests a 
bulge membership. The PN PK 167-0.1 is projected close to 
the OC New Cluster 1, and the similar distances are con- 
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